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(1) .  Administration  information 

Grant  Number:  FA9550-10-1-0533 
CRYOGENIC  PELTIER  COOLING 

Ohio  State  University/Princeton  University/Massachusetts  Institute  of 
Technology/Boston  College/Michigan  State  University/California  Institute  of 
Technology/University  of  Houston 

Project  Director:  Joseph  P.  Heremans,  Ohio  State  University  (Ohio  Eminent  Scholar, 
Mechanical  Engineering,  Physics) 

Princeton  PI:  Robert  Cava  (Chemistry) 

Massachusetts  Institute  of  Technology  Pis:  Mildred  S.  Dresselhaus  (Physics),  Gang  Chen 
(Mechanical  Engineering) 

Boston  College  PI:  Cyril  Opeil  (Physics) 

Michigan  State  University  PI:  Donald  T.  Morelli  (Materials  and  Chemical  Engineering) 
California  Institute  of  Technology:  G.  Jeffrey  Snyder  (Materials  Science) 

University  of  Houston  PI:  Zhifeng  Ren  (Physics) 

(2) .  Program  objective/scientific  approach 

This  program  is  aimed  at  developing  the  fundamental  science  and  technology  to 
provide  solid-state  cooling  in  the  10  K  to  150  K  temperature  range.  The  ultimate 
temperature  difference  that  each  stage  of  a  Peltier  cooler  can  achieve,  and  the 
thermodynamic  efficiency  of  a  Peltier  cooler,  are  a  function  of  essentially  one  material 
parameter,  the  thermoelectric  figure  of  merit  zT=T  S2cj/k ;  where  T  is  the  average 
temperature  of  operation,  S  is  the  thermopower  or  Seebeck  coefficient,  and  cr  and  k  are 
the  electrical  and  thermal  conductivity  respectively.  The  goal  of  the  program  is  to 
conduct  research  into  new  materials  and  new  physical  phenomena  that  will  lead  to  an 
enhancement  of  zT. 

The  difficulty  is  that  the  materials  parameters  S,  cr  and  k  are  mutually  counter- 
indicated,  i.e.  any  change  made  to  increase  one  almost  always  results  in  a  decrease  of 
another:  increasing  the  electron  concentration  to  increase  cr  decreases  S  by  a  fundamental 
relation,  called  the  Pisarenko  relation.  Increasing  the  degree  of  disorder  in  the  solid  to 
decrease  /c  almost  always  decreases  cr.  The  problem  involves  the  study  of  transport 
properties  of  both  the  electrons,  which  dominate  S  and  cr,  and  of  the  phonons,  which 
dominate  k. 

The  scientific  approaches  used  in  this  MURI  are  categorized  into  two  classes: 

(1)  Increase  S  for  a  given  cr  by  band  structure  engineering:  indeed,  it  is  known  that 
a  sharp  local  increase  at  a  specific  energy  of  the  electronic  density  of  states  (DOS)  leads 
to  an  optimal  product  S2cr  [G.D.  Mahan  and  J.O.  Sofo,  Proc.  Natl  Acad.  Sci.  93,  7436 
(1996)].  This  is  done  specifically  by  using  (la)  resonant  dopants  in  semiconductors; 
(lb)  semiconductors  based  on  transition  metals,  where  the  d-levels  of  select  atoms  create 
the  desired  DOS;  (lc)  heavy  Fermion  metals,  where  the /-levels  of  rare  earths  create  the 
desired  DOS;  and  (Id)  size  quantization. 

(2)  Decrease  k  for  a  given  cr  by  nanostructuring,  thus  creating  defects  of  one 
particular  length  scale  so  as  to  scatter  phonons  but  not  electrons.  This  is  done  (2a)  by 
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developing  heat- treatments  during  the  casting  of  ingots  of  TE  semiconductors;  (2b)  by 
mechanical  nano  structuring  using  ball-milling  and  Spark  Plasma  Sintering. 

The  materials  that  were  the  focus  of  the  initial  efforts  (but  these  have  evolved  over  the 
duration  of  the  program),  were  selected  for  having  the  best  zT  at  the  onset  of  the  program: 

(1)  Bii-xSbx  alloys,  very  narrow-gap  semiconductors.  The  rc-type  materials  had 
the  best  zT  <0.5  at  the  onset  of  the  program  [B.  Lenoir  et  al.,  J.  Phys.  Chem.  Solids  59, 
129  (1998).]  and  the  p- type  materials  reached  zT  <  0.08  at  the  onset  of  the  program.  [Y.S. 
Hor  and  R.J.  Cava,  J.  Alloys  Compd.  479,  368  (2009).] 

(2)  Tetradymite  semiconductors,  of  the  general  formula  (Bii-xSbx)2(Tei-ySey)3, 
which  are  used  in  commercial  Peltier  coolers,  and  in  the  six-stage  Peltier  coolers  that 
reach  the  lowest  cryogenic  temperatures  to  date  (about  180  K). 

(3)  Heavy  Fermion  metals  [G.D.  Mahan  and  J.O.  Sofo,  Proc.  Natl.  Acad.  Sci.  93, 
7436  (1996).] 

(4)  Semiconductors  with  transition  metal  atoms  as  principal  constituents,  initially 
FeSb2  [A.  Bentien  et  al.,  EPL  80,  17008  (2007).] 

(3) ,  The  most  recent  program  review  (date,  location,  and  government 
participants) 

Title:  “Cryogenic  Peltier  Cooling  4th  Year  Review  Meeting” 

Date:  3  May  2015  -  5  May  2015 

Location:  Kirtland  Air  Force  Base,  Albuquerque,  NM 

Government  attendees:  Dr.  Harold  Weinstock  AFOSR/RSE;  Dr.  Charles  Stein 

AFOSR/Kirtland  AFB  (ret.);  Col.  William  Cooley,  Chief  of  the  Space  Vehicles 

Directorate,  AFRL;  Douglas  Dudis,  WPAFB;  Keith  Avery,  Kirtland  AFB;  William 

Byrne,  Kirtland  AFB. 

MURI  team  members  in  attendance:  Joseph  P.  Heremans  (OSU);  Robert  Cava 
(Princeton);  Mildred  S.  Dresselhaus  and  Gang  Chen  (MIT);  Cyril  Opeil  (BC);  Zhifeng 
Ren  (UH);  Donald  T.  Morelli  (MSU);  G.  Jeffrey  Snyder  (Caltech);  numerous  post¬ 
doctoral  fellows  and  graduate  students  supported  by  the  program. 

(4) ,  Financial  execution 

The  period  of  performance  and  funding  is  as  follows: 

Initial  Award  Increment  #1  (2  months)  9/14/10-11/14/10  $271,257 
Increment  #2  (12  months)  11/15/10-11/14/11  $1,500,000 
Increment  #3  (12  months)  11/15/11-11/14/12  $1,500,000 
Increment  #4  (10  months)  11/15/12  -  9/14/13  $1,  228,743 
Increment  #5  (2  months)  9/14/13  -  11/14/13  $294,537 
Increment  #6  (12  months)  11/15/13  -  11/14/14  $1,500,000 
Increment  #7  (10  months)  11/15/14-9/14/15  $1,205,463 

TOTAL  FUNDING  TO  DATE:  $7,500,000 
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In  this  MURI  program  as  each  increment  of  funds  is  added,  the  funds  are  split  to  be 
approximately  OSU  25%,  Princeton  12.5%,  Michigan  State  12.5%,  Cal  Tech  12.5%, 
Boston  College/University  of  Houston  18.75%  (split  28.5%/71.5%),  and  MIT  18.75%. 
Expenditure  rates  are  consistent  with  the  rate  of  incremental  funding. 


(5),  Accomplishments  to  date: 

5a  Transport  theory 
5a,  1  Electron  transport 

Heavy  Fermion  systems  are  believed  to  be  good  candidates  for  thermoelectric 
cooling  applications  at  low  temperature  because  of  the  sharp  features  of  density  of  states 
near  the  Fermi  level,  which  gives  rise  to  good  Seebeck  coefficients.  In  combination  with 
nano  structuring  process,  which  drastically  reduces  the  lattice  thermal  conductivity,  such 
systems  hold  promise  of  possessing  high  ZT  values.  To  understand  the  effects  of 
nano  structuring  on  electron  transport  in  heavy  Fermion  systems,  thermoelectric 
properties  in  a  Kondo  insulator  nano  wire  were  investigated  by  the  Chen  group  at  MIT. 
The  bulk  electronic  transport  are  simulated  within  the  framework  of  Dynamic  Mean  Field 
Theory  (DMFT),  showing  very  small  electron  mean  free  paths,  which  indicates  the 
electron  properties  can  be  preserved  in  nanostructures,  while  the  phonons  with  longer 
mean  free  paths  can  be  strongly  scattered  by  the  interfaces  due  to  classical  size  effects. 
Both  effects  contribute  to  high  ZT  values  up  to  1.5  even  at  low  temperatures.  The  results 
have  been  published  by  Nano  Letters. 

We  are  investigating  the  potential  of  using  modulation  doping  and  anti-resonance 
scattering  in  core-shell  nanoparticles  to  improve  the  electron  performance.  Initial 
simulation  shows  that  this  is  a  promising  direction. 

Experimentally,  a  characterization  system  was  built  to  study  electron  transport  in 
thermoelectric  material  below  room  temperature.  The  setup  is  capable  of  measuring 
electrical  conductivity,  Seebeck,  Hall,  and  isothermal  Nernst  coefficients  from  4K  to 
350K.  The  measured  transport  coefficients  can  be  analyzed  to  find  energy  dependence  of 
electron  scattering,  density-of- states  effective  mass,  Fermi  level  and  electron’s  mobility. 
To  validate  the  system,  nano-grained  Bi2Te2.7Seo.3  sample  was  tested  and  the  measured 
electrical  conductivity  and  Seebeck  coefficients  match  well  with  ones  measured  by  a 
commercial  characterization  system  (ULYAC  ZEM3)  at  300K  with  discrepancy  less  than 
5%.  This  experimental  setup  will  be  applied  to  nano-grained  BiSb  alloys. 

5a,2  Lattice  thermal  conductivity  from  first  principles 

We  have  applied  first  principles  calculation  to  simulate  phonon  thermal 
conductivity  and  mean  free  path  in  materials  of  interests  to  our  program,  including  Bi, 
Sb,  and  FeSb2. 

Bi-Sb  alloy  is  one  of  the  best  thermoelectric  materials  for  refrigeration 
applications.  To  control  coupled  transport  of  electrons  and  phonons  and  to  enhance  zT 
further,  it  is  essential  to  understand  their  transport  in  the  material.  However,  in  bismuth 
and  Bi-Sb  alloy,  phonon  contribution  to  the  total  thermal  conductivity  is  not  well  known 
since  electrons  and  phonons  have  comparable  contributions  to  total  thermal  conductivity. 
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To  quantify  phonon  thermal  conductivity  and  understand  phonon  transport  in  bismuth, 
we  calculated  lattice  dynamics  in  bismuth  using  first  principles  and  the  Boltzmann 
equation  which  was  very  accurate  in  the  case  of  other  thermoelectric  materials  such  as  Si- 
Ge  alloy.  The  calculated  lattice  thermal  conductivity  matches  very  well  with  the  previous 
experimental  work  in  low  temperature  range.  Also,  phonon  mean  free  path,  which  is 
important  for  designing  nanostructure  sizes  to  enhance  phonon  scattering,  was 
successfully  calculated.  From  the  calculation,  we  found  three  interesting  features  of 
phonon  transport  in  bismuth.  First,  the  9th  nearest  neighbors  in  the  structure  cannot  be 
ignored  even  though  the  interatomic  distance  is  long.  This  peculiar  feature  originates 
from  the  rocksalt-like  crystal  structure  of  bismuth  and  large  electronic  polarizability,  and 
significantly  affects  both  harmonic  and  anharmonic  interactions  between  atoms.  Second, 
the  previous  work  about  lattice  thermal  conductivity  of  bismuth  estimated  that  lattice 
thermal  conductivity  is  much  smaller  than  electronic  thermal  conductivity  near  room 
temperature.  However,  according  to  our  calculation,  lattice  thermal  conductivity  is 
comparable  to  the  electronic  thermal  conductivity  near  room  temperature.  Therefore,  the 
nanostructure  approach  to  reduce  thermal  conductivity  by  enhanced  phonon  scattering 
may  be  still  effective  for  bismuth.  Last,  bismuth  is  well-known  for  its  anisotropic  electron 
transport  properties,  but  phonon  transport  in  bismuth  is  almost  isotropic.  This  may  be  due 
to  small  distortion  of  its  crystal  structure  from  exact  rocksalt  structure. 

FeSb2  has  been  intensively  studied  recently  as  a  promising  candidate  for  low 
temperature  thermoelectric  applications  for  its  extremely  large  Seebeck  coefficient 
around  10K,  while  the  thermoelectric  efficiency  is  limited  by  its  high  thermal 
conductivity.  We  have  been  studying  the  thermal  transport  properties  of  FeSb2,  including 
the  thermal  conductivity  and  especially  the  phonon  mean  free  path,  from  first-principles 
simulations,  aiming  at  providing  a  guideline  to  reduce  the  thermal  conductivity  of  FeSb2 
by  nano  structuring.  We  extracted  the  second  and  the  third  order  force  constants  between 
atoms  in  the  lattice  from  the  force-displacement  data  calculated  using  density  functional 
theory  (DFT).  The  information  of  the  second  order  force  constants  enabled  the 
calculation  of  the  phonon  dispersion,  which  was  in  good  agreement  with  that  calculated 
using  alternative  density  functional  perturbation  theory  (DFPT).  Anharmonic  properties 
of  FeSb2  can  be  deduced  from  the  third  order  force  constants,  based  on  which  we 
obtained  the  temperature-dependent  thermal  expansion  coefficient  that  matches  well  with 
the  experimental  data.  To  validate  the  third  order  force  constants,  we  calculated  the  mode 
Grlineisen  parameters  from  the  extracted  force  constants  and  compared  with  those 
calculated  using  DFPT.  Good  agreement  was  observed  and  the  calculated  thermal 
conductivity  was  found  in  decent  agreement  with  experimental  data  when  the  sample 
boundary  scattering  was  also  taken  into  account  at  low  temperatures.  The  phonon  mean 
free  paths  were  also  extracted,  which  indicated  there  is  enough  room  for  nanostructuring 
to  reduce  phonon  thermal  conductivity. 
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5b  Experimental  investigations  of  Heavy  Fermion  metals 

The  Morelli  group  at  MSU  has  been  interested  in  exploring  thermoelectric  effects 
below  room  temperature  in  alloys  comprised  of  elements  containing  f-shell  electrons.  It 
is  well  known  [G.D.  Mahan  and  J.O.  Sofo,  Proc.  Natl.  Acad.  Sci.  93,  7436  (1996).]  that  a 
delta-function-type  electron  transport  distribution  function  can  maximize  the 
thermoelectric  properties  by  providing  simultaneously  large  electrical  conductivity  o  and 
large  Seebeck  coefficient  a.  While  such  a  distribution  is  not  rigorously  found  in  nature, 
it  can  be  approximately  achieved  in  alloys  comprised  of  elements  with  f-shell  valence 
electrons;  these  electrons  form  a  narrow  band  in  the  density  of  states  near  the  Fermi  level. 
We  have  chosen  to  focus  initially  on  two  families  of  compounds:  i)  those  based  on 
CePd3;  and  ii)  those  based  on  YbAh.  The  primary  reasons  for  this  choice  are  firstly,  that 
these  two  base  compounds  possess  the  largest  power  factors  (a2a)  of  known  compounds 
in  the  cryogenic  temperature  regime,  and  secondly  they  offer  both  p-type  and  n-type 
behavior,  respectively,  a  necessary  requirement  for  the  development  of  a  Peltier  cooler. 

A  large  part  of  our  work  up  to  this  point  has  focused  on  identifying  and 
implementing  mechanisms  for  reducing  the  lattice  thermal  conductivity  by  compositional 
and  structural  modification,  while  at  the  same  time  maintaining  large  power  factor.  We 
have  explored  three  separate  approaches  for  increasing  phonon  scattering  in  these 
materials: 

•  Solid  solution  formation.  This  includes  Cei-xScxPd3,  CePd3  xPtx,  Ybi-xScxAh,  and 
Ybi-x(RE)xAb,  where  RE  =  Lu  or  Er; 

•  “ Resonant  ”  scattering  by  filler  atoms.  This  includes  CePd3Mx,  with  M  =  Al,  Ga, 
Si,  and  Ge; 

•  Nano  structure -induced  phonon  scattering  in  CePd3. 

5b.l  CePcb-based  materials 
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1.  The  CePcb  crystal  structure  (isostructural 
with  AuCu3)  is  shown  in  Figure  1 .  The  total  thermal 
conductivity  is  a  sum  of  electronic  and  lattice  (phonon) 
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Figure  2.  Thermoelectric  figure  of 
merit  in  the  CdPch-xPtx  solid  solution. 

The  value  of  >  0.  is  the  world  record 
for  a  p-type  material  at  100  K. 

that  the  lattice  thermal  conductivity  is  practically  unchanged,  a  behavior  suggesting  that 
the  Ce  site  plays  a  negligible  role  in  phonon  propagation.  Substituting  Pd  for  Pt,  on  the 
other  hand,  gives  rise  to  the  classic  reduction  in  lattice  thermal  conductivity.  At  the  same 
time,  the  Seebeck  coefficient  is  hardly  changed  at  all  with  substitution,  behavior  we 
attribute  to  the  close  ion  size  match  between  Pd  and  Pt,  allowing  for  the  preservation  of 
the  intermediate  valent  state  that  gives  rise  to  large  Seebeck  coefficient.  As  a  result 
(Figure  2),  the  figure  of  merit  ZT  =  (a2a)T/K,  is  increased  for  the  solid  solution  alloys 
and  reaches  values  rivaling  the  best  p-type  materials  below  room  temperature.  [S.R. 
Boona  and  D.T.  Morelli,  Appl.  Phys.  Lett.  101,  101909  (2012).] 

A  second  and  distinct  approach  [S.R.  Boona  and  D.T.  Morelli,  J.  Appl  Phys.  112, 
063709  (2012).]  investigates  structural  modification  by  forming  alloys  of  the  form 
CePdsMx  (0<  x  <  1).  Here  the  M  atom  occupies  a  vacant  body-centered  void  position  in 
the  AuCu3  lattice  (Figure  1).  Such  “filling”  of  large  voids  in  structures  is  known  to 
strongly  affect  the  physical  properties,  and  especially  the  thermoelectric  properties,  the 
skutterudites  being  a  case  in  point.  We  find,  again,  a  strong  reduction  in  thermal 
conductivity  that  saturates  with  both  filling  fraction  and  filler  mass.  Unfortunately,  the 
filled  structures  invariably  exhibit  depressed  Seebeck  coefficient,  and  the  figure  of  merit 
is  enhanced  only  slightly  at  100K. 


contributions;  the 
latter  can  be 
modified  by  solid 
solution 

formation,  as 
differences 
between  the  mass 
and  size  of 
substituted  atoms 
can  give  rise  to 
increased  phonon 
scattering.  Since  ScPd3  is  isostructural  with  CePd3, 
and  Sc  and  Ce  have  a  large  mass  difference,  one 
might  expect  solid  solutions  of  Cei-xScxPd3  would 
exhibit  lattice  thermal  conductivity  reduction. 
Instead,  on  samples  with  x  ranging  from  zero  to 
0.25,  we  find  [S.R.  Boona  and  D.T.  Morelli, 
Journal  of  Electronic  Materials  41,  1199  (2012).] 


Figure  1.  The  YbAbMx  structure. 
Large  red  circles:  Yb;  large  open 
circles:  Al;  small  filled  circle:  M 
ions. 


2.  Finally,  we  have  investigated  whether  nanostructuring  of  CePd3  can  lead  to 
enhanced  thermoelectric  properties,  primarily  through  thermal  conductivity  reduction  via 
interface  scattering  of  phonons.  [S.R.  Boona  and  D.T.  Morelli  (in  preparation).]  Again, 
this  approach  results  in  very  effective  reduction  in  thermal  conductivity.  However,  the 
delicate  intermediate  valence  state  giving  rise  to  high  Seebeck  coefficient  is  perturbed  in 
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nanostructured  materials,  likely  due  to  an  increasing  fraction  of  trivalent  Ce  at 
nanocrystal  surfaces  and  grain  boundaries. 

3.  An  important  outcome  of  the  CePd3-based  work  is  the  development  of  some 
basic  strategies  for  reducing  thermal  conductivity  while  maintaining  high  Seebeck 
coefficient.  [S.R.  Boona  and  D.T.  Morelli,  Journal  of  Electronic  Materials  DOI 
10. 1007/sl  1664-012-2328-7  (2012).]  For  instance,  Seebeck  coefficients  and  power 
factors  can  be  enhanced,  though  modestly,  using  chemical  pressure,  temperature,  or 
doping;  our  continuing  and  future  efforts  are  focused  on  combining  these  effects  with 
some  of  the  thermal  conductivity  reduction  mechanisms  described  above  in  materials 
systems  similar  to  CePd3. 


5b.2  YbAh-and  YbAh-based  materials 

Efforts  in  the  YbAL  compound  family 
have  primarily  focused  on  the  influence  of 
substitutions  on  the  Yb  site.  Solute  ions  included 
Sc,  Lu,  and  Er.  Representative  results  are 
described  here  for  the  Ybi-xScxAl3  solid  solution. 

We  devised  a  novel  method  for  synthesis  of  these 
compounds  that  involves  a  combination  of  arc¬ 
melting,  annealing,  etching,  and  powder 
processing  that  allows  for  the  fabrication  of  this 
alloy  family  for  the  first  time.  [G.J.  Lehr  and  D.T. 

Morelli,  Intermetallics  32,  225  (2013).] 

Substitution  of  Sc  results  in  a  drastic  reduction  in 
electrical  resistivity  and  Seebeck  coefficient;  the 
resistivity  of  pure  ScAb  is  on  the  order  of  that  of 
pure  Al,  and  the  Seebeck  coefficient  is  near  zero. 

Interestingly,  however,  for  intermediate  values  of  x,  the  Seebeck  coefficient  drops  much 
less  rapidly  with  Sc  concentration,  and  as  a  result,  at  intermediate  compositions,  the 
power  factor  is  enhanced  over  pure  YbAL.  At  the  same  time,  the  increase  in  electronic 
thermal  conductivity  is  counterbalanced  by  the  expected  reduction  in  lattice  thermal 
conductivity  due  to  alloy  scattering.  Taken  together,  this  leads  to  an  increase  of  the 

figure  of  merit  by  some  30%  (Figure  3). 
Similar  results  are  obtained  using  Lu 
and  Er.  [G.J.  Lehr  and  D.T.  Morelli, 
Journal  of  Electronic  Materials  DOI 
10. 1007/sl  1664-012-2401-2  (2012).] 

Finally,  because  use  of  this  material  in  a 
Peltier  module  will  require  mechanical 
property  data,  we  have  undertaken  a 
study  [R.D.  Schmidt,  E.  D.  Case,  G.  J. 
Lehr,  and  D.  T.  Morelli,  Intermetallics 
35,  15  (2013).]  of  the  elastic  modulus, 
shear  modulus,  and  Poisson’s  ratio  of 
Figure  4.  Thermoelectric  power  factor  in  the  Ybi-  these  compounds  and  studied  the  effects 
xScxAb  solid  solution 
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Figure  3.  Thermoelectric  figure  of 
merit  in  the  Ybi-xScxAl3  solid  solution 
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of  porosity  in  powder  processed  samples  on  these  mechanical  parameters. 

In  parallel  with  the  investigations  on  YbAL,  we  have  also  explored  thermoelectric 
effects  in  the  Ybi-xScxAh  family.  [G.J.  Lehr  and  D.T.  Morelli  (in  preparation).]  This 
alloy  is  of  interest  because,  unlike  the  case  of  YbAL,  substitution  of  Sc  in  YbAh 
increases  the  lattice  constant  and  drives  the  Yb  ions  to  a  more  intermediate  valent  state 
and  thus  has  the  potential  for  increasing  the  Seebeck  coefficient.  Our  preliminary 
investigations  on  this  alloy  family  show  that  indeed  this  is  the  case.  As  can  be  seen  in 
Figure  4,  both  the  Seebeck  coefficient  (not  shown)  and  the  power  factor  increase 
dramatically  as  the  Sc  concentration  is  increased  up  to  90%,  before  falling  slightly; 
additionally,  the  temperature  of  the  peak  value  in  the  power  factor  shifts  to  lower 
temperature  with  increasing  Sc.  Both  of  these  observations  are  consistent  with  a 
narrowing  of  the  Kondo  resonance  peak  in  the  f-electron  density  of  states  of  this  system. 
It  is  interesting  to  note  that  large  power  factor  does  not  require  large  numbers  of  Yb  ions 
-  rather,  only  that  they  be  spatially  separated  sufficiently  to  maximize  the  Kondo 
resonance.  This  result  is  interesting  and  potentially  important  as  it  gives  us  a  means  of 
tuning  both  the  magnitude  of  the  power  factor  and  the  temperature  at  which  the 
maximum  occurs,  which  can  be  made  to  coincide  with  the  required  temperature  of 
operation  of  a  Peltier  cooler.  While  the  ZT  of  this  system  does  not  exceed  0.1,  it  is  hoped 
that  this  concept  can  be  applied  to  other  intermediate  valence  systems  with  larger 
potential  figures  of  merit. 


5c  Bii-xSbx  alloys 

Based  on  their  bulk  properties,  bismuth  antimony  (BiSb)  alloys  with  low  Sb 
concentration  are  expected  to  be  the  materials  showing  the  best  currently  known 
thermoelectric  (TE)  properties  at  temperatures  below  77  K,  the  temperature  of  liquid 
nitrogen.  The  work  on  these  alloys  falls  in  three  categories: 

5c.  1  Bulk  and  single  crystal  alloys  doped  with  resonant  impurities,  which  is  further 
broken  down  into  work  on  n-type  and  work  on  p-type  alloys 
5c. 2  Quantum  confined  structures 
5c.3  Bulk  nanostructured  composites 

5c.l  Bulk  and  single  crystal  alloys  doped  with  resonant  impurities 

5c.l.l  n-type  single  crystalline  Bii-xSbx  alloys  (OSU): 

First  principles  calculations  on  various  impurities  in  Bi  were  done  by  our 
collaborator,  Bartlomiej  Wiendlocha  at  AGH  University  of  Science  and  Technology 
(Krakow,  Poland).  Among  those  impurities,  Potassium  (K)  was  found  to  be  a  promising 
resonant  impurity  in  Bi  while  Lithium  (Li)  was  predicted  to  be  an  n-type  dopant  if  it  goes 
into  interstitial  sites.  Figure  1(a)  shows  that  K  introduces  a  sharp  peak  in  density  of  states 
(DOS)  of  Bi  near  £>  which  resembles  the  known  resonant  T1  peak  in  PbTe.  [C.M. 
Jaworski,  B.  Wiendlocha,  V.  Jovovic  and  J.P.  Heremans,  Energy  and  Environ.  Sci .,  4(10), 
4155  (2011).]  From  this  result,  it  is  expected  that  if  K  impurity  level  is  hybridized  into 
the  extended  state  of  host  Bi,  it  will  increase  DOS  of  Bi  where  its  sharp  peak  is  located. 
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When  Ef  hits  the  distorted  DOS,  the  resonant  effect  will  appear  in  transport  properties. 
Fig.  1(b)  indicates  that  Li  tends  to  increase  Ef  when  it  takes  the  interstitial  sites  in  Bi. 


Figure  1.  First  principles  calculations  for  Density  of  States  of  Bi  with  (a)  K  impurity  and  (b)  Li  impurity 
in  interstitial  sites 

The  band  structure  of  Bii-xSbx  alloys  is  very  similar  to  that  of  Bi  except  that  Bii- 
xSbx  with  0.07<x<0.22  has  a  band  gap  while  Bi  has  the  band  overlap  between  the 
conduction  band  at  the  L-point  and  the  valence  band  at  the  T-point.  Therefore,  the 
behavior  of  K  and  Li  in  Bii-xSbx  alloys  is  anticipated  to  be  similar  to  that  in  Bi.  The 
calculation  of  the  DOS  of  Bii-xSbx  alloys  with  K  and  Li  impurities  is  underway.  Fig.  2 
shows  the  experimental  data  on  thermoelectric  (TE)  properties  of  K  doped  Bii-xSbx 
samples  compared  to  the  undoped  sample.  It  is  obvious  that  K  introduces  an  anomaly  in 
S,  exhibiting  an  additional  hump  at  20K  and  much  larger  S  at  lower  temperatures.  For 
Bi93Sb7  +  0.1%  K  sample,  the  power  factor  is  increased  by  almost  2  times  compared  to 
the  undoped  sample  and  as  a  result,  zT  is  enhanced  about  50%  from  100K  up  to  300K, 
reaching  the  maximum  ~0.7  at  120K  which  is  a  record  value  up  to  date. 


T  (K) 


T  (K) 


Figure  2.  Seebeck  coefficient  ( S ),  Electrical  resistivity  (p),  Power  factor,  and  zT  of  Big^Sbio.s  (Orange), 
Bi93Sb7  +  0.1  %K  (Green),  and  BLoSbio  +  1%K  (Dark  green)  samples.  The  values  ZT>0.5  set  the  world 
record  for  n-type  materials  in  this  temperature  range. 

S  and  p  of  the  undoped  Bi8sSbi2  and  1%  Li  doped  Bi8sSbi2  single  crystal  samples 
are  presented  in  Fig.  3.  [C.M.  Orovets,  A.M.  Chamoire,  H.  Jin,  B.  Wiendlocha  and  J.P. 
Heremans,  Journal  of  Electronic  Materials  41,  1648  (2012).]  The  Li  doped  sample 
shows  metallic  behavior  compared  to  the  undoped  sample.  These  experimental  results 
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confirm  that  Li  actually  goes  into  the  interstitial  sites  in  BissStm  so  it  acts  as  a  strong  n- 
type  dopant. 


(Red) 
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5c.l.2  p-type  single  crystalline  Bii-xSbx  alloys: 

There  are  two  main  results  obtained  so  far  for  p-type  single  crystalline  Bii-xSbx 
alloys.  Firstly,  Indium  (In)  was  found  to  be  a  new  acceptor  in  Bi  and  Bii-xSbx  alloys 
which  is  surprising  since  In  is  the  same  trivalent  element  as  Bi.  Fig.  4  shows  that 
Shubnikov  de-Haas  (SdH)  measurements  reveal  that  In  increases  the  hole  concentration 
in  Bi. 

(b) 


B(T) 

Figure  4.  (a)  Shubnikov  de-Haas  (SdH)  oscillations  in  the  electrical  resistivity  for  pure  Bi,  Bi+0.1%  In, 
and  Bi+0.5%In  samples.  The  oscillations  are  periodic  in  1/H  where  H  is  the  magnetic  field  intensity,  (b) 
Analysis  results  from  SdH  oscillations  from  (a)  where  Af  is  the  cross  sectional  area  of  the  Fermi  surface,  Ef 
the  Fermi  level,  and  Nh  the  hole  carrier  concentration.  It  is  obvious  that  In  introduces  additional  holes  in 
Bi. 


Samples 

[Ad/H)]1 

(T) 

Af 

(m-2) 

Ef 

(meV) 

Nh  (cm'3) 

Bi 

6.352 

6.06  xlO16 

11.5 

2.93  x  1017 

Bi  +  0.1%  In 

7.665 

7.31  xlO16 

13.9 

3.88  x  1017 

Bi  +  0.5%  In 

12.09 

1.15  x  1017 

21.9 

7.69  x  1017 

0  2  4  6  8 


Secondly,  it  was  shown  [H.  Jin,  C.M.  Jaworski  and  J.P.  Heremans,  Appl.  Phys. 
Lett.  101,  053904  (2012).]  that  doping  multiple  valence  bands  can  enhance  the  TE 
performance  of  p-type  Bii-xSbx  alloys  (Figure  5). 


1 1 .6<Sb  <26.5  30<Sb<50  70<Sb<90 


T  (K) 


Figure  5.  (a)  Pisarenko’s  plot  of  thermopower  versus  hole  concentration  at  T=80K  for  all  measured 


samples  in  this  study.  Solid  line  and  dashed  line  are  calculated  for  the  valence  bands  at  the  H-points  and 
the  valence  band  at  the  T-point,  respectively.  Points  are  experimental  data.  Carrier  density  of  each  sample 
was  obtained  from  Hall  measurement,  (b),  (c),  (d)  Schematic  band  diagrams  for  the  three  different  Sb 
concentration  regimes  shown,  (e)  Samples  with  highest  zT  are  presented.  Those  samples  correspond  to  the 
case  (b)  and  show  thermopower  above  the  calculated  Pisarenko’s  line. 
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5c.2  Quantum  confined  Bii-xSbx  alloys  structures  at  MIT 

Since  further  enhancements  of  TE  properties  comes  from  quantum  confinement 
effects,  further  TE  enhancement  is  expected  in  thin  film  BiSb  alloys,  which  as  films  also 
show  the  best  three  dimensional  packing  for  nano-structured  materials.  However,  the 
transport  properties  of  BiSb  alloys  are  very  anisotropic.  For  these  reasons,  we  have  spent 
the  past  2.5  years  calculating  various  aspects  of  the  electronic  structure  of  BiSb  thin  films 
as  a  function  of  crystal  orientation,  thickness  and  temperature  up  to  77K.  These 
calculations  are  now  being  extended  to  finding  and  optimizing  the  pertinent 
thermoelectric  properties.  With  the  start  of  the  final  two  years  of  the  program,  we  have 
added  a  postdoctoral  researcher  to  grow  appropriate  samples  and  to  make  thermoelectric 
measurements  as  a  function  of  Sb  concentration  in  the  BiSb  alloys,  their  crystal 
orientation,  film  thickness  and  temperature.  The  following  text  gives  details  on  the 
ongoing  theoretical  and  experimental  programs  and  their  coupling. 


5c.2.1  Theoretical  Work  on  quantum  confined  Bii-xSbx 

The  Dresselhaus  group  has  theoretically  studied  the  electronic  band  structure 
properties  of  BiSb  thin  films  as  a  function  of  film  growth  orientation,  film  thickness,  and 
stoichiometry  for  the  cryogenic 


temperature  range.  We  have  developed 
a  general  model  for  the  low¬ 
dimensional  narrow-band-gap  material 
systems,  where  only  the  lowest 
conduction  band  and  the  highest 
valence  band  are  important  for  the 
transport  properties.  The  model  for  a 
BiSb  thin  film  system  is  an  iterative- 
two-dimensional-two-band  model. 
Firstly,  the  two-band  model  describes 
the  narrow-band-gap  and  the  non¬ 
parabolic  dispersion  relation  in  bulk  Bi 
and  bulk  BiSb  alloys  with  remarkable 
accuracy.  Because  of  the  special 
properties  of  Bi  and  BiSb  materials,  the 
de  Broglie  wavelength  of  electrons  in 
these  materials  is  ultra-long,  so  that 
quantum  size  effects  become  obvious 
even  when  the  film  thickness  is  as  large 
as  100  nm.  Thus,  it  is  difficult  to  carry 
out  first  principles  calculations,  such  as 
density  functional  theory  calculations, 
to  model  the  Bi  and  BiSb  thin  films,  due 
to  the  heavy  computation  caused  by 
both  the  nano-  and  alloy-effects.  To 
solve  this  problem,  we  developed  an 
iterative  approach  to  use  the  two-band 
model  to  study  the  two-dimensional 
(2D)  system.  The  narrow-band-gap  and 


Trigonal 
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100  200 
4  (nm) 

Fig.  1.  (a)  How  the  L(1)-point  single-Dirac -cone  is 
formed.  The  figure  shows  the  band  structure  of  a 
Bio.96Sbo.04  thin  film  grown  normal  to  the  bisectrix  axis 
and  how  the  band  structure  changes  over  different  film 
thicknesses.  The  green  curves  show  the  lowest 
conduction  band  (upper  one)  and  the  highest  valence 
band  (lower  one)  at  the  L(2)  and  L(3)  points  in  the 
Brillouin  zone.  The  blue  curves  are  for  the  L(1)  point. 
The  dashed  red  curve  is  the  highest  valence  band  at  the 
T  point.  The  L(1)-point  band  gap  remains  less  than 
~1  meV  until  the  film  thickness  lz  <  80  nm.  The  L(2)  and 
L(3)  points  have  the  same  band  gap,  which  is  largely 
opened  up.  Thus,  an  anisotropic  single-Dirac-cone  is 
formed  at  the  L(1)  point  when  the  film  thickness  lz  is 
large  enough  to  retain  the  L(1)-point  mini-gap  at  values 
less  than  ~1  meV.  (b)  How  to  control  the  electronic 
anisotropy  of  the  Dirac -cone  systems  of  Bii-xSbx  thin 
films.  A  general  scheme  of  how  the  Fermi  velocity  (vf) 
anisotropy  coefficient  y  =  lvmaxl/lvminlfor  the  L(1)-point 
single  Dirac  cone  changes  with  film  growth  orientation. 
The  value  of  y  for  a  specific  film  growth  orientation  is 
shown  by  the  radius  and  the  color  on  the  left  scale,  y 
can  be  as  large  as  ~14,  and  as  small  as  ~2.  According  to 
our  calculations,  y  is  sensitive  neither  to  the  film 
thickness  /,  nor  to  the  Sb  composition  x  in  the  range  of 
0  <  x  <  0.15. 
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the  non-parabolic  dispersion  relations 
are  obtained  when  the  iterations  come  to 
convergence.  This  iterative-two- 
dimensional  two-band  model  may  be 
extended  to  a  more  general  iterative 
low-dimensional-N-band  model,  if 
needed. 

From  our  simulations,  we 
theoretically  point  out  the  possibility  of 
constructing  different  types  of  Dirac- 
cone  materials  based  on  BiSb  thin  films, 
as  seen  in  Fig.  1.  We  found  that  we  can 
engineer  the  band  structure  of  BiSb  as  a 
function  of  the  film  stoichiometry, 
growth  orientation,  and  thickness.  We 
predict  that  it  is  possible  to  construct 
Dirac  cones  with  different  degrees  of 
anisotropy  and  topology  numbers,  such 
as  single-Dirac-cone  materials,  bi-Dirac- 
cone  materials  and  tri-Dirac-cone 
materials  as  in  Fig.  2a-c.  Furthermore, 
we  also  predict  how  to  construct  quasi- 
Dirac-cone  materials  with  different 
magnitude  mini-band-gaps,  and  how  to 
construct  semi-Dirac- cones,  where  the 
electrons  are  linearly  dispersed  in  one 
direction  and  parabolically  dispersed  in 
the  other  direction,  as  described  in  Fig. 
2d. 


I(2)  Lm  L(3> 


. . .  r  i  i  i  if  'i  i  ii 

-5  0  5-5  0  5-5  0  5 


k  (nnr1) 


Fig.  2.  An  illustration  of  (a)  single-,  (b)  bi-  and  (c)  tri- 
Dirac-cone  Bii-xSbx  thin  films  grown  respectively 
along  the  (a)  bisectrix,  (b)  binary  and  (c)  trigonal  axes. 
The  illustration  is  based  on  the  example  of  Bii-xSbx  thin 
films  with  film  thickness  /  =  100  nm,  x  =  0.04, 
P  =  1  atm  and  T  <11  K,  under  which  the  L  points  of 
bulk  Bii-xSbx  have  a  zero-gap. (d)  Example  of  a  semi- 
Dirac  cone  in  the  Bii-xSbx  thin  film  system  (x  =  0.10 
and  /  =  100  nm).  We  see  that  around  the  L(1)  point,  the 
fermions  are  linearly  dispersed  along  the  vmax  direction, 
and  parabolically  dispersed  along  the  Vmin  direction. 


To  guide  thin  film  crystal 

growers,  we  have  drawn  out  the  electronic  phase  diagrams,  in  Fig.  3,  as  a  function  of  film 
growth  orientation,  film  thickness,  and  stoichiometry,  based  on  our  theoretical 
calculations,  which  include  semi-metal  (SM)  phases,  direct  semi-conductor  (DSC)  phases 
and  indirect  semi-conductor  phases  (ISC).  Finally,  we  have  calculated  the  overall  band 
gap  and  band  overlap,  as  a  function  of  film  growth  orientation,  film  thickness  and 
stoichiometry  in  Fig.  3d-f.  Furthermore,  to  control  the  properties  of  these  cryogenic 
thermoelectric  materials,  we  calculated  the  carrier  concentration  as  a  function  of  film 
orientation,  film  thickness,  stoichiometry,  and  Fermi  level  in  Fig.  4,  to  study  their 
different  transport  properties.  Future  theoretical  work  will  be  focused  on  modifying  our 
model  to  account  for  the  measured  results  and  for  calculating  thermoelectric  properties  in 
detail. 
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Fig.  3.  (a-c)  Phase  diagrams  of  Bii-xSbx  thin  films  for  different  growth  orientations  as  a  function  of  film 
thickness  and  stoichiometry.  The  regions  of  semi-metal  (SM),  indirect-semiconductor  (ISC),  and  direct- 
semiconductor  (DSC)  phases  are  shown  explicitly,  (d-f)  The  diagrams  of  overall  narrow-band-gap  and 
band  overlap  of  Bii-xSbx  thin  films  for  different  growth  orientations  as  a  function  of  film  thickness  and 
stoichiometry.  The  magnitude  of  the  band  gap  for  a  semiconductor  phase  is  positive,  while  the 

magnitude  of  the  band  overlap  for  a  semi-metal  phase  is  negative. 


5c.2.2  Experimental  Work  on 
quantum  confined  Bii-xSbx 

Experimental  verification  of 
various  theoretically  predicted 
Dirac  cones  in  BiSb  hinges  on 
controlling  the  thin  film  crystal 
orientation.  Since  the  films  are  Bi 
rich,  we  focus  first  on  studying 
the  growth  dynamics  of  Bi. 

Previous  techniques  of 
fabricating  Bi  and  BiSb  thin  films  made  samples  with  the  trigonal  axis  perpendicular  to 
the  substrate,  by  lattice  matching  Bi  to  a  supporting  substrate  (e.g.  Mica  or  BaF2).  The 
key  to  growing  high  quality  thin  films  is  the  choice  of  substrate  and  growth  conditions 
(e.g.  deposition  technique,  temperature,  and  pressure). 

Thin  Bi  films  were  deposited  on  Mica  using  both  thermal  evaporation  (TE)  and 
molecular  beam  epitaxy  (MBE).  Figure  5a  and  5b  shows  atomic  force  microscopy  images 
of  both  types  of  films  (~20  nm  thick).  The  TE  deposited  film  has  grain  boundaries  that 
are  shaped  as  triangles,  which  are  characteristic  of  films  grown  along  the  trigonal 
direction.  The  MBE  grown  films,  however,  produce  a  much  more  diverse  array  of  grain 
boundary  shapes.  This  indicates  that  there  is  a  mixture  of  crystallographic  orientations. 
We  use  MBE  to  study  the  impact  of  the  substrate  on  the  growth  because  MBE  provides  a 
more  diverse  set  of  grains.  We  grew  Bi  films  on  Si  (111),  Si  (001),  and  sapphire 
substrates.  The  sapphire  substrate  stood  out  as  having  larger  and  more  triangular  grains 
than  either  Si  substrate  (Fig.  5c  and  5d).  To  investigate  the  differences,  we  took  x-ray 
diffraction  (XRD)  spectra  (see  Fig.  5e).  From  these  spectra  we  identify  two  different 
preferred  crystal  orientations.  One  is  the  trigonal  (001)  (coordinates  for  Bi  are  given 
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Fig.  4.  Carrier  concentration  of  Bio.96Sbo.04  vs.  film  thickness 
lz  and  Fermi  level  Ef  are  shown  at  (a)  77  K  and  (b)  4.2  K. 
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using  rhombohedral  indices,  as  shown 
in  Fig.  5)  and  the  other  is  along  the 
(012)  direction.  From  the  relative 
heights  of  the  (001)  and  (012)  peaks, 
we  see  that  sapphire  produces  more 
trigonal  oriented  grains  when 
compared  with  Bi  films  grown  on  Si. 
This  result  is  expected  since  sapphire 
is  better  lattice  matched  to  Bi  than  Si. 
Still,  the  results  demonstrate  that  the 
choice  of  substrate  based  on  the  lattice 
mismatch  can  lead  to  growth  of 
different  crystal  orientations. 

5c.2.3  Future  directions  of  the  work 
on  Bi  films  and  quantum-confined 
structures 

Since  the  substrate  plays  a 
critical  role  in  defining  the  crystal 
orientation  of  the  film,  we  will  fine 
tune  the  substrate  lattice  parameters. 
BaF2  has  been  shown  to  be  an 
excellent  match  to  trigonal  oriented  Bi 
films.  To  change  the  lattice  constant 
of  BaF2  itself,  we  can  grow  a  different 
material  and  use  it  as  a  buffer  layer 
between  BaF2  and  our  supporting 
substrate.  The  gradually  applied  strain 
from  the  buffer  layer  will  slightly 
change  the  lattice  constant  of  the 
BaF2. 


Fig.  5.  Atomic  force  microscopy  (AFM)  images  of 
-20  nm  thick  Bi  films,  (a)  Bi  film  deposited  by  thermal 
evaporation  (TE)  produces  triangular  grains 
(characteristic  of  the  trigonal  axis)  on  the  order  of 
-100  nm.  (b)  Bi  film  grown  by  molecular  beam  epitaxy 
(MBE)  produces  grains  of  various  shapes  also  on  the 
order  of  -100  nm.  A  few  grains  are  outlined  to  guide 
the  eye.  AFM  images  of  MBE  grown  Bi  thin  films  on  a 
(c)  Si  (111)  and  (d)  sapphire  substrate.  The  Bi  grains 
on  sapphire  are  larger  and  consist  of  more  triangles,  (e) 
X-ray  diffraction  spectrum  of  Bi  film  grown  on 
sapphire,  Si  (111),  and  Si  (001).  The  two  primary 
crystal  orientations  of  the  Bi  in  a  rhombohedral  basis 
( ai,a2,c )  are  the  trigonal  face  (001),  and  (012).  The 
inset  shows  a  close-up  of  the  (024)  peak,  indicating 
that  Bi  grown  on  sapphire  exhibits  the  smallest  peak 
here,  since  its  grains  are  primarily  oriented  along  the 
trigonal  axis. 


Once  samples  are  produced,  we  will  use  transport  measurements  and 
spectroscopic  techniques  to  search  for  grains  with  quasi-Dirac  cones.  Gated  structures 
will  be  fabricated,  which  will  allow  us  to  sweep  the  Fermi  level  through  the  electronic 
bands.  All  thermo-electric  relevant  properties  will  also  be  measured.  Thermal 
conductivity  can  be  measured  using  optical  reflectance  techniques.  The  Seebeck 
coefficient  can  be  measured  by  using  micro-fabricated  platforms  to  simultaneously  bias 
the  device  while  monitoring  temperature  and  heat  flow.  Finally,  all  results  will  be  used  to 
calibrate  our  theoretical  models.  This  will  allow  us  to  theoretically  optimize  film 
thicknesses,  operating  temperatures,  and  doping  concentrations.  Theoretical  modeling 
will  be  coordinated  with  experimental  findings. 
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5c.3  Nanostructuring  of  Bii-xSbx  at  Caltech 

The  reduction  of  thermal  conductivity,  specifically  lattice  thermal  conductivity, 
has  been  a  technique  used  to  increase  the  figure  of  merit  in  thermoelectric  materials  for 
over  twenty  years.  The  reduction  is  achieved  by  scattering  of  phonons  at  the  boundary 
between  two  different  structures.  However,  this  reduction  is  only  beneficial  if  it  is 
achieved  while  maintaining  a  power  factor  of  approximately  the  same  value  or  larger 
than  if  the  structures  did  not  exist.  Here,  the  Caltech  group  focused  on  ways  to 
incorporate  a  second  phase  into  alloys  of  Bi-Sb  in  order  to  reduce  the  lattice  thermal 
conductivity  without  significantly  hindering  the  electron  transport  in  the  system. 

Prior  studies  of  Bi-Sb  alloys  exploited  the  use  of  grain  size  reduction  to  increase 
the  overall  figure  of  merit.  The  aim  of  this  work  was  to  use  a  metallurgical  approach  to 
improve  the  thermoelectric  properties  of  Bi9oSbio.  Concurrently,  using  bulk  processing 
and  consolidation  techniques,  this  should  act  to  improve  the  mechanical  stability  of  these 
easily  cleaved  layered  compounds  as  well. 

Attempts  using  SnTe,  Ag,  and  As  have  been  made,  with  the  most  success 
involving  precipitation  of  As  samples.  Typically,  single  crystal  measurements  are  done 
both  in  the  trigonal  direction  and  perpendicular  to  in  order  to  maximize  zT,  however  we 
have  focused  on  a  subject  not  often  studied,  that  is,  the  properties  of  polycrystalline  Bi- 
Sb.  The  previous  studies  of  polycrystalline  Bi-Sb  focusing  on  grain  size  reduction  have 
never  shown  great  promise  because  of  the  fact  that  the  electron  mobility  is  severely 
limited.  However,  the  goal  here  is  to  maintain  a  large  grain  size,  with  a  small  enough 
structure  size  to  adequately  scatter  phonons  without  debilitating  the  electron  mobility. 

In  the  past,  two  types  of  consolidation  have  been  used.  The  first  involves  cold 
pressing  and  then  sintering  the  material  with  the  second  involving  hot  extrusion  of  the 
cold  pressed  material.  The  density  of  often  far  below  that  of  the  theoretical  density  for 
samples  that  are  cold  pressed  and  sintered,  and  often  times  not  reported  in  the  literature. 
However,  samples  that  were  extruded  after  cold  pressing  had  densities  which  approached 
that  of  the  theoretical  value.  These  samples  were  shown  to  have  preferential  orientation 
relative  to  the  extrusion  direction,  resulting  in  texturing  and  directionally  dependent 
properties.  We  have  been  working  on  utilizing  our  rapid  induction  heated  hot  press  to 
achieve  similar  results  to  that  of  the  extruded  samples,  but  with  the  benefit  of  having  little 
to  not  texturing. 

Experimental  Results 


Figure  1.  Scanning  electron  microscope  (SEM)  image  in  backscatter  mode  of  Bi-SnTe  sample  melted  and 
annealed 
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Initial  attempts  to  nanostructure  elemental  bismuth  with  SnTe  were  unsuccessful 
due  to  the  low  temperatures  necessary  for  the  precipitation  step.  Typically,  materials  are 
melted,  homogenized  and  then  precipitated  in  three  separate  isothermal  steps.  However, 
since  Bi  has  such  a  low  melting  temperature  and  the  reaction  kinetics  are  slow  enough  at 
the  temperatures  in  which  it  was  believed  to  be  homogeneous  that  this  attempt  was 
unsuccessful. 


Figure  2.  Electron  microscope  images  of  BkoSbio  with  Ag:  (a)  melt,  (b)  2  at.  %,  and  (c)  4  at.  % 
microstructure 


As  can  be  seen  in  Figures  3  and  4,  both  attempts  at  nano  structuring  Bi9oSbio  were 
successful  in  creating  second  phase  with  the  potential  to  effectively  scatter  phonons. 
Figure  5a  shows  the  melt  micro  structure  for  the  Ag  added  Bi-Sb  sample.  After  annealing 
for  3  days  at  250°C  both  the  2  at.  %  (Fig.  6b)  and  4  at.  %  (Fig.  7c)  Ag  samples  had 
similar  microstructure.  However,  because  of  the  doping  effect  of  Ag  in  Bi-Sb  alloys 
further  measurements  were  not  conducted. 


Figure  8.  Electron  microscope  images  of  BigoSbio  with  As:  (a)  melt,  (b)  2  at.  %,  and  (c)  4  at.  % 
microstructure 


The  As  structured  samples,  however,  had  a  marked  difference  between  the  2  at.  % 
(Fig.  9b)  and  4  at.  %  (Fig.  10c)  microstructure.  The  large  precipitates  seen  in  Figure  11a 
are  formed  during  melting  step,  however  there  is  still  a  significant  amount  of  As 
dissolved  in  the  matrix  that  then  phase  segregates  upon  annealing  at  250°C  for  three  days. 
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Figure  12.  Electron  microscope  images  of  (a)  consolidated  BigoSbio  and  (b)  with  4  at.  %  As 
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There  is  an  issue  with  precipitate  coarsening  during  the  rapid  consolidation  step  as 
seen  in  Fig.  13.  The  consolidation  parameters  work  well  with  regard  to  densification  as 
the  sample  hot  pressed  without  arsenic  showed  little  to  no  porosity,  with  a  density  at  or 
above  98%  of  the  theoretical  (Fig.  14a).  However,  the  4  at.  %  sample  that  was  hot 
pressed  showed  a  large  amount  of  material  coarsening  (Fig.  15b).  The  fine  micro  structure 
seen  in  the  inset  of  Fig  16c  had  coarsened  and  all  of  the  fine  microstructure  was  no  longer 
present. 


Figure  17.  Transport  measurements  for  pure  BboSbio  perpendicular  and  parallel  to  the  pressing  direction  to 
establish  the  extent  of  texturing. 


The  resulting  low  temperature  transport  measurements  revealed  that  the  sample 
anisotropy  seems  to  be  minimal  (Fig.  18).  The  property  most  problematic  with  regard  to 
anisotropy  is  the  thermal  conductivity  (Fig.  19c)  as  there  is  a  drastic  difference  in  the 
room  temperature  thermal  conductivity.  Based  on  the  low  temperature  measurements 
conducted,  this  anisotropy  is  significantly  reduced  indicating  little  to  no  preferential 
orientation  during  hot  pressing  with  both  samples  having  comparable  zT  (Fig.  20d). 


Temperature  (K) 


Figure  21.  Low  temperature  transport  measurements  for  pure  BkoSbio,  2  at.  %  and  4  at.  %  arsenic 
structured  samples. 


The  first  round  of  low  temperature  transport  measurements  revealed  that  the 
arsenic  structures  do  indeed  reduce  the  thermal  conductivity  as  expected  (Fig.  22c). 
However,  the  results  of  the  Seebeck  and  resistivity  measurements  (Fig.  23a-b)  are 
difficult  to  explain  solely  from  the  transport  measurements  as  each  sample  had  an 
equivalent  or  near  equivalent  carrier  concentration.  It  is  believed  that  arsenic  could  be 
increasing  the  value  of  the  band  gap,  but  more  work  is  required  to  determine  if  this  is  the 
case. 


Conclusion 

Successful  micro  structure  formation  has  been  attained  in  alloys  of  Bi-Sb. 
The  current  method  for  micro  structure  generation  has  proven  successful  at  thermal 
conductivity  reduction  for  samples  structured  in  the  range  of  2-4  at.  %  arsenic.  However, 
it  is  believe  that  these  results  can  be  improved  by  eliminating  the  precipitate  coarsening 

18 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


upon  consolidation.  This  can  potentially  be  remedied  by  annealing  the  samples  for  less 
time  prior  to  consolidation.  This  is  something  that  will  be  examined  in  the  near  future. 

Once  it  has  been  established  that  the  coarsening  can  be  reduced  during 
consolidation,  continued  coordination  with  the  Ohio  State  group  for  the  purpose  of 
conducting  low  temperature  transport  measurements  will  be  necessary  to  assess  the 
extent  to  which  the  zT  can  be  improved  for  the  optimized  samples. 

Also  to  be  examined  is  the  possibility  of  band  gap  widening  due  to  the  addition  of 
arsenic.  This  will  be  examined  with  respect  to  the  lattice  parameter  change  as 
characterized  by  x-ray  diffraction. 


5d  Mechanical  nanostruturing,  particularly  of  FeSb2 

The  MURI  project  activities  at  Boston  College  are  a  collaborative  effort  between 
the  Opeil  Group,  which  specializes  in  thermodynamic  material  characterization  at 
cryogenic  temperatures  and  the  Ren  Group  responsible  for  sample  materials  preparation. 
The  main  thrust  of  our  research  combines  the  technique  of  nanostructured  material 
preparation  and  the  intriguing  report  of  a  large  Seebeck  coefficient  (-45,000  pV/K)  in 
single  crystal  FeSb2.  [Sun  et  al.,  Dalton  Trans.  39,  1012  (2010).]  Our  research  efforts 
have  focused  on  improving  the  thermoelectric  performance  ( ZT)  of  FeSb2  using  the 
techniques  of  chemical  doping,  the  integration  of  Cu  and  Ag  nanoparticles  [H.  Zhao,  et 
al.,  Nanotechnology  23,  505402  (2012).],  nano  structuring  [H.  Zhao  et  al.,  Appl.  Phys. 
Lett.  99,  163101  (2011).],  and  stoichiometric  adjustment.  [M.  Pokharel  et  al.  MRS 
Proceedings  1456,  1505  (2012).]  This  research  has  led  to  significant  improvements  in 
thermoelectric  efficiency,  ZT,  primarily  through  a  significant  reduction  of  thermal 
conductivity  resulting  from  increased  phonon  scattering  and  a  decrease  in  resistivity 
resulting  from  the  introduction  of  Ag  nanoparticles  during  the  nanocomposite  synthesis. 
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Figure  1.  Comparison  of  the  temperature  dependent  thermoelectric  efficiency,  ZT,  of  nanocomposite  FeSb2 
with  different  grain  sizes  and  FeSb2  single  crystal  data  from  Sun,  et  al.  J.  Phys.  Conf.  Ser.  150,  012049 
(2009).  The  optimized  nanocomposite  sample  represents  an  increase  of  160%  over  the  single  crystal 
performance,  for  further  details  see  H.  Zhao  et  al.,  Appl.  Phys.  Lett.  99,  163101  (2011). 
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Figure  2.  Temperature  dependent  thermoelectric  efficiency,  ZT,  comparison  for  an  increasing  amount  of 
Ag  nanoparticles  in  the  FeSt)2  nanocomposite  matrix.  The  optimized  ZT  occurs  for  FeSb2Agm  at  m  =  0.11 
due  to  the  metal/semiconductor  interface  engineering  approach  stemming  from  the  similar  work  functions 
of  FeSb2  and  Ag.  This  approach  reduced  both  the  thermal  conductivity  and  resistivity.  For  details  see,  H. 
Zhao  et  al.,  Nanotechnology  23,  505402  (2012). 


In  the  following,  the  activities  and  accomplishments  for  the  first  two  and  a  half 
years  are  summarized. 

Note  on  Professor  Zhifeng  Ren’s  move  to  the  University  of  Houston 

As  part  of  the  MURI  program  on  Cryogenic  Cooling  using  thermoelectric 
materials  led  by  Prof.  Joseph  Heremans  at  Ohio  State  University,  Ren’s  group  has  carried 
out  work  on  materials  synthesis  and  characterizations  at  Boston  College  (BC)  before  the 
end  of  2012.  Starting  from  January  1  2013,  his  group  moved  from  BC  to  University  of 
Houston  (UH).  In  the  past  two  and  a  half  years,  his  group  has  first  doubled  the  ZT  of  the 
single  crystals  reported  in  literature  in  the  first  year  by  creating  nanostructured  FeSb2 
using  ball  milling  and  hot  pressing  technique  his  group  developed  over  the  last  a  few 
years,  then  doubled  the  doubled  ZT  in  the  second  year  by  adding  nanoparticles  into  the 
nanostructured  FeSb2.  At  UH,  his  group  will  continue  to  carry  out  the  proposed  work 
with  enhanced  intensity  due  to  better  facilities  his  group  is  setting  up  at  UH. 

5d.l  Preparation  and  characterization  of  nanostructured  FeSb2 

Nanostructured  FeSb2  samples  are  made  through  ingot  synthesis  by  melting,  ball 
milling,  and  hot  pressing  processes.  By  optimizing  the  processing  parameters, 
nanostructured  FeSb2  materials  showed  the  lowest  thermal  conductivity  of  less  than  1 
W/mK  that  is  much  lower  than  the  peak  thermal  conductivity  of  500  W/mK  of  the  single 
crystals,  however  the  Seebeck  coefficient  also  decreased  drastically  to  about  125  pV/K, 
only  a  tiny  fraction  of  the  single  crystals’  42,000  pV/K.  Nevertheless,  a  thermoelectric 
figure-of-merit  (ZT)  of  0.013  was  achieved,  which  is  more  than  double  of  the  state-of- 
the-art  value  of  the  single  crystals.  Interestingly  the  ZT  curve  is  also  dramatically 
broadened,  which  is  very  beneficial  to  the  applications  at  a  broader  temperature  range.  In 
addition,  doping  in  this  material  has  also  been  studied. 
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5d.2  Preparation  of  nanostructured  Bi-Sb  alloy  materials 

The  alloy  was  made  by  melting  method.  The  nanostructured  materials  were 
prepared  by  ball  milling  and  hot  pressing.  In  addition,  melt  spinning  was  also  used  to 
synthesize  Bi-Sb  ribbons  with  the  size  of  micrometers.  In  the  process,  we  have  added 
potassium  (K)  into  the  materials.  Interesting  results  were  found  in  the  K-doped  Bi-Sb 
polycrystalline  materials.  An  enhancement  of  about  20%  at  300  K  on  power  factor  from 
both  the  Seebeck  coefficient  and  electrical  conductivity  was  observed. 

5d.3  Fabrication  of  nanostructured  FeSb2-xAgx  with  nanoinclusions  Agi-vSbv 

Based  on  the  result  of  enhanced  thermoelectric  figure-of-merit  ( ZT)  of  0.013  on 
nanostructured  FeSb2  we  achieved  in  2011,  we  tried  to  further  improve  the  ZT  of 
nanostructured  FeSb2  through  doping.  Out  of  the  many  elements  we  have  tried  so  far,  we 
found  that  silver  (Ag)  provided  some  enhancement  to  the  nanostructured  FeSb2.  By 
replacing  some  Sb  with  Ag,  we  synthesized  nanostructured  FeSb2-xAgx  with  10  %  (molar 
ratio)  nanoinclusions  Agi-ySby  with  significant  ZT  enhancement  comparing  to  the 
nanostructured  pure  FeSb2.  The  experimental  data  show  that  both  thermal  conductivity 
and  electrical  resistivity  of  the  nanocomposites  were  significantly  reduced  in  the  lower 
temperature  regime  comparing  with  pure  FeSb2  owing  to  the  interfaces  between  FeSb2- 
xAgx  and  Agi-ySby  phases,  as  well  as  the  similar  work  functions  between  the  two  phases. 
Overall,  significant  improvement  in  ZT  (from  0.013  to  0.019)  was  achieved  for  the 
optimized  nanocomposite  composition  FeSb1.975Ag0.025/Ag0.77Sb0.23  in  which  Ago.77Sbo.23 
is  about  10%  by  molar  ratio.  The  results  by  this  approach  clearly  demonstrated  that 
metal/semiconductor  interface  played  an  important  role  and  confirmed  the  potential  of 
strongly  correlated  material  systems  as  promising  thermoelectric  materials. 

5d.4  ZT  enhancement  of  nanostructured  FeSb2  by  adding  Cu  nanoparticles 

After  the  work  done  on  Ag  addition  into  FeSb2,  we  realized  that  Cu  may  be  even 
better  than  Ag  since  Cu  does  not  react  with  FeSb2  whereas  Ag  does.  The  isolated  Cu 
nanoparticles  can  act  as  modulation  dopant  to  provide  the  needed  electrons  in  increasing 
the  electrical  conductivity  without  degrading  the  Seebeck  coefficient.  Sure  enough, 
addition  of  Cu  nanoparticles  have  made  us  to  achieve  nanocomposite  FeSb2Cuo.o45  with 
ZTs  of  0.026,  such  a  ZT  value  doubled  our  last  year’s  record  of  0.013,  fulfilling  the 
promise  of  ZT  doubling  every  year  since  we  started  this  MURI  program  in  2010. 

5d.5  Investigation  of  other  thermoelectric  materials  working  at  cryogenic 
temperatures 

Aside  from  FeSb2,  we  are  also  trying  to  find  new  materials,  such  as  CrSb2,  which 
has  been  reported  to  having  Seebeck  coefficient  of  5000  |iV/K  below  10  K.  At  the  same 
time,  we  are  also  trying  to  synthesize  another  p-type  material  CeCu6. 

5d.6  Analysis  of  Phonon  Drag  in  FeSb2  Nanocomposite 

In  an  attempt  to  understand  the  physics  of  the  large  Seebeck  coefficient  in  FeSb2, 
reported  by  Sun  et  al.  [Sun  et  al.,  Dalton  Trans ,  39,  1012  (2010).]  we  studied  the  thermal 
conductivity  (k)  and  Seebeck  coefficient  (S)  of  nanocomposites  of  varying  grain  size. 
The  grain- size  dependence  of  these  properties  and  the  correlations  between  the  properties 
themselves  reveals  that  phonon-drag  contributes  significantly  to  the  large  S -value  of 
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coarse-grained  samples.  The  severe  degradation  of  S-peak  for  nanocomposites  is  a 
consequence  of  the  suppression  of  the  phonon-drag  effect  due  to  increased  boundary 
scattering. 


T(  K) 


Figure  3.  Temperature  dependence  of 
the  Seebeck  coefficient  for  the  four 
different  grain  size  samples.  Fitting  to 
the  power  law  was  done  for  all  the 
samples  below  50  K.  Inset:  Grain  size 
dependence  of  the  peak  value  (Smax) 
and  the  peak  position  (Tmax)  of  the 
Seebeck  coefficient.  Smax  for  all 
samples  is  significantly  small  than  that 
of  a  single  crystal  at  low-T,  whereas  at 
high  T  the  values  are  comparable.  As 
shown  in  the  inset,  the  Smax  decreases 
with  decreasing  grain  size.  This  is 
expected  for  the  phonon-drag  effect, 
because  the  non-electronic  scattering 
reduces  the  phonon  mean  free  path, 
which  reduces  the  phonon-drag 
contribution. 


Our  nanocomposite  data  and  analysis  of  phonon  drag  do  not  support  other 
mechanisms  that  might  explain  the  large  Seebeck  effects  reported  by  Bentien  et  al.  [A. 
Bentien  et  al.,  Euro.  Phys.  Lett.  80,  17008  (2007).]  Although  one  cannot  preclude  the 
presence  of  electron-electron  correlation  effects,  their  role  in  this  phenomenon  may  be  a 
minor  one.  The  recent  analysis  of  electron  correlations  using  a  hybrid  functional 
approach  of  Becke  [A.D.  Becke,  J.  Chem.  Phys.  98,  1372  (1993).]  and  Hegin’s  GW 
functional  approach  [L.  Hedin,  Phys.  Rev.  139,  A796  (1965).]  by  Tomczak  et  al.  [J.M. 
Tomczak,  Phys.  Rev.  B  82,  085104  (2010).]  suggest  that  the  high  thermopower  in  FeSb2 
should  not  be  understood  in  the  context  of  local  correlations,  but  rather  by  utilizing  vertex 
corrections  to  the  transport  coefficients.  Such  vertex  corrections  describe  the  phonon- 
drag  effect.  The  phonon-drag  effects  in  FeSb2  are  similar  to  those  described  in  p-type  Ge 
by  Geballe  and  Hull.  [T.H.  Geballe  and  G.W.  Hull,  Phys.  Rev.  94,  1134  (1954).]  In  a 
similar  vein  the  study  of  magnetoresistance  and  Hall  effect  by  Takahashi  et  al.  [H. 
Takahashi  et  al.,  Phys.  Rev.  B  84,  205215  (2011).]  concludes  that  the  large  Seebeck 
coefficient  in  FeSb2  is  unlikely  to  originate  from  electron-electron  correlations  because 
they  have  an  insignificant  effect  on  the  Seebeck  coefficient  in  the  low -temperature 
insulating  regime.  Our  data  on  FeSb2  nanocomposite  supports  their  conjecture  that  the 
phonon-drag  effect  plays  an  essential  role  in  enhancing  the  Seebeck  coefficient  in  the  low 
temperature  regime,  as  shown  in  other  semiconductor  materials  e.g.,  InSb  [H.P.R. 
Frederikse  and  E.V.  Mielczarek,  Phys.  Rev.  99,  1889  (1955).]  and  weakly  P  doped  Si.  [L. 
Weber  and  E.  Gmelin,  Appl.  Phys.  A  53,  136  (1991).] 

5d.7  Kapitza  Resistance  and  Thermal  Conductivity  in  FeSb2  Nanocomposite 

In  further  analyzing  the  strong  grain- size  dependent  thermal  conductivity  of  the 
FeSb2  nanocomposites,  we  used  a  simple  macroscopic  model  (effective  medium 
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approximation,  EMA)  [D.A.G.  Bruggeman,  Ann.  Phys.  (Leipzig)  24,  636  (1935);  R. 
Landaure,  J.  Appl.  Phys.  23,  779  (1952).]  based  on  Kapitza  resistance  calculation  using 
equation,  (T,d)  =  — »  where  tq  is  the  intrinsic  (or  bulk)  thermal  conductivity,  Rk  is 

i+  §  1 
a 

the  Kapitza  resistance  and  d  is  the  grain  size.  Here  Ki  is  grain-size  independent  but  is 
temperature  dependent  while  Lk  —  RkKt  and  Gk  =  —  are  called  the  Kapitza  length  and 

Rk 

Kapitza  conductance,  respectively.  In  an  effective  medium  approach,  the  total  thermal 
conductivity  of  a  polycrystal  is  the  sum  of  the  intra-grain  and  the  inter-grain 
contributions.  The  inter-grain  component  arises  due  to  the  interfacial  resistance,  also 
known  as  the  Kapitza  resistance,  to  the  thermal  transport.  In  the  presence  of  thermal 
gradient,  the  Kapitza  resistance  results  in  a  temperature  discontinuity  across  the  interface, 
which  was  first  reported  by  Kapitza  [P.L.  Kapitza,  J.  Phys.  (Moscow)  4,  181  (1941)]. 
The  overall  effect  of  the  Kapitza  resistance  is  to  reduce  the  thermal  conductivity  and  the 
effect  is  more  pronounced  at  smaller  grain  sizes.  Our  analysis  showed  that  the  Kapitza 
length  found  for  nanocomposite  FeSb2  at  10  K  was  390  nm  with  a  Kapitza  resistance  of 
the  order  of  10"7  m2KW4,  which  confirmed  the  importance  of  interfacial  thermal 
resistance  across  the  grain  boundaries  in  determining  the  thermal  properties  of 
nanocomposite  FeSb2  sample.  Our  thermal  conductivity  data  shows  the  Kapitza 
resistance  was  the  predominating  factor  in  determining  the  bulk  thermal  properties  below 
50  K. 


T(K) 


Figure  4.  Temperature  dependence 
of  the  Kapitza  resistance  (Rk)  and 
Kapitza  conductance  (Gk)  obtained 
as  best  fit  to  equation  given  above. 
Rk  and  Gk  values  are  represented  by 
left  and  right  Y-axes,  respectively. 

^  The  open  symbols  with  dashed  lines 
represent  the  temperature 
dependence  of  the  bulk  thermal 
of  resistance  (d/xi)  as  indicated  by  the 
left  Y-axis  for  the  five  samples. 
Stars:  S-5000,  Diamonds:  S-300, 
down-pointed  triangles:  sample  S- 
100,  up-pointed  triangles:  sample  S- 
30,  open-circles:  sample  S-20,  filled 
circles:  Gk  (T),  filled  sq.:  Rk  (T). 


5e  The  Thomson  cooler  derived  from  the  Compatibility  Factor  Approach  to  the 

analysis  of  Thermoelectric  Devices  (Caltech) 

Traditional  thermoelectric  Peltier  coolers  exhibit  a  cooling  limit  which  is 
primarily  determined  by  the  figure  of  merit,  zT.  Rather  than  a  fundamental 
thermodynamic  limit,  this  bound  can  be  traced  to  the  difficulty  of  maintaining 
thermoelectric  compatibility.  Self-compatibility  locally  maximizes  the  cooler’s 
coefficient  of  performance  for  a  given  zT  and  can  be  achieved  by  adjusting  the  relative 
ratio  of  the  thermoelectric  transport  properties  that  make  up  zT.  We  investigated  the 
theoretical  performance  of  thermoelectric  coolers  that  maintain  self-compatibility  across 
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the  device.  We  find  that  such  a  device  behaves  very  differently  from  a  Peltier  cooler,  and 
we  term  self- compatible  coolers  “Thomson  coolers”  when  the  Fourier  heat  divergence  is 
dominated  by  the  Thomson,  as  opposed  to  the  Joule,  term.  A  Thomson  cooler  requires  an 
exponentially  rising  Seebeck  coefficient  with  increasing  temperature,  while  traditional 
Peltier  coolers,  such  as  those  used  commercially,  have  comparatively  minimal  change  in 
Seebeck  coefficient  with  temperature.  When  reasonable  material  property  bounds  are 
placed  on  the  thermoelectric  leg,  the  Thomson  cooler  is  predicted  to  achieve 
approximately  twice  the  maximum  temperature  drop  of  a  traditional  Peltier  cooler  with 
equivalent  figure  of  merit  (zT). We  anticipate  that  the  development  of  Thomson  coolers 
will  ultimately  lead  to  solid-state  cooling  to  cryogenic  temperatures.  In  order  to  make 
such  a  cooler  new  kinds  of  high  Seebeck  and  low  Seebeck  materials  are  required.  This 
opens  up  new  challenges  and  opportunities  for  thermoelectric  material  development. 

5f  Tetradymites 

At  the  outset,  the  bulk  of  our  research  of  the  Cava  group  at  Princeton  focused  on 
improving  the  thermoelectric  properties  of  state-of-the-art  bismuth  telluride  selenide 
materials,  which  we  attempted  through  several  different  approaches.  However,  by  the 
present  point  in  time  we  have  shifted  our  focus  towards  a  few  other  classes  of  materials 
as  well,  as  discussed  below.  With  regard  to  bismuth  telluride  selenide,  we  primarily 
concentrated  on  alloys  that  are  selenium-rich  (i.e.,  Bi2Te2Se  and  Bi2TeSe2)  compared  to 
the  tellurium-rich  state-of-the-art  materials  that  have  already  been  optimized  through 
countless  other  studies  [see  e.g.,  D.M.  Rowe,  CRC  Handbook  of  Thermoelectrics,  CRC 
Press,  1995].  Because  relatively  Se-rich  materials  have  received  less  attention  in  the 
literature,  their  defect  chemistry  is  less  understood  and  consequently  less  control  over 
their  properties  has  been  achieved.  Thus,  our  studies  on  these  materials  have  sought  to 
understand  and  control  the  defect  chemistry  and  correspondingly  the  carrier 
concentrations  in  these  materials,  which  directly  affect  the  thermoelectric  properties. 

The  work  presented  in  [Shuang  Jia  et  al.,  Phys  Rev  B  84,  235206  (2011)]  reports 
the  characterization  of  the  material  BhTe2Se,  which  was  synthesized  by  two  different 
techniques  -  i.e.,  the  modified  Bridgman  and  Bridgman-Stockbarger  (B-S)  methods  -  and 
with  different  amounts  of  Te-excess/Se-deficiency.  These  crystals  yielded  different  sets 
of  electronic  properties  based  on  different  concentrations  of  defects;  in  particular,  it  was 
noteworthy  that  crystals  of  low  carrier  concentrations  (down  to  ~5  x  1014  cm'3  at  4  K) 
were  achieved  and  could  be  grown  consistently.  Incidentally,  this  finding  was  of  interest 
to  the  topological  insulator  community  (which  is  also  intently  studying  Bi-Te-Se 
compounds  at  present).  Yet,  the  ability  to  achieve  such  a  large  extent  of  compensation  in 
this  material  is  also  noteworthy  from  the  perspective  of  our  thermoelectrics  research,  and 
so  we  have  pursued  further  investigations  on  this  material.  Specifically,  we  have  spent 
the  past  year  probing  the  homogeneity  of  B-S  grown  as  well  as  Sn-substituted  BhTe2Se 
crystals,  and  have  collected  a  significant  amount  of  thermoelectric  data  that  have  been 
drafted  into  a  manuscript  [M.K.  Fuccillo,  Shuang  Jia,  M.E.  Charles,  and  R.J.  Cava, 
Thermoelectric  properties  of  BETe2Se  compensated  by  native  defects  and  Sn  doping] 
recently  submitted  to  the  Journal  of  Electronic  Materials.  An  additional  study  [Shuang  Jia 
et  al.,  Phys  Rev  B  86,  165119  (2012)]  that  we  performed  on  BETe2Se  involved 
collaboration  with  the  Princeton  University  Department  of  Physics  to  use  scanning 
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tunneling  microscopy  to  “see”  and  identify  the  defects  in  crystals  of  slight  Bi-excess/Te- 
deficiency. 

Another  direction  we  set  out  on  while  studying  the  bismuth  tellurides  involved  the 
material  Bi2Te1.6S1.4j  which  exhibits  rich  defect  chemistry  and  electronic  properties 
similar  to  those  of  Bi2Te2Se  because  of  the  stronger  electronegativity  of  S  atoms  (like  Se) 
compared  to  Te.  Although  we  supported  characterization  of  the  material  properties  of  this 
compound  form  a  thermoelectrics  perspective,  Bi2Tei.6Si.4  turned  out  to  be  more  of 
interest  as  a  topological  insulator  [Huiwen  Ji  &  al.,  Phys  Rev  B  85  201103(R)  (2012)] 
than  a  thermoelectric  material. 

Yet  another  avenue  we  have  pursued  on  Bi-Te-Se  compounds  was  to  consider  the 
material  Bi2TeSe2.  The  study  outlined  in  [M.K.  Fuccillo  &  al.,  Solid  State 
Communications  152  1208  (2012)]  shows  how  we  used  Sb  substitution  for  Bi  to  hole- 
dope  and  compensate  a  strongly  n-type  material  to  the  point  of  an  n-p  crossover  and  then 
beyond  it  into  the  heavily-doped  p-type  side.  We  studied  the  effects  of  this  substitution 
on  the  thermoelectric  properties  and  found  that  the  thermoelectric  figure  of  merit  over  the 
composition  range  studied  was  highest  near  the  crossover. 

As  is  evident  from  the  above,  much  of  our  work  on  the  bismuth  tellurides  has  supported 
topological  insulator  research.  The  influence  of  this  work  on  that  field  has  recently  been 
discussed  in  a  review  article  submitted  to  the  Journal  of  Materials  Chemistry  [R.J.  Cava 
&  al.  (Submitted  to  the  Journal  of  Materials  Chemistry)]. 

The  Cava  group  has  recently  begun  to  shift  focus  towards  studying  FeSb2-like  and 
MnTe2-like  compounds  as  described  in  the  MURI  proposal  and  section  5d  above.  They 
have  determined  several  different  possible  crystal  growth  methods  to  obtain  pure  and 
doped  crystals  of  many  of  these  materials,  and  a  number  of  crystals  have  been 
successfully  obtained.  Electronic  transport  measurements  and  Investigations  into  the 
defect  or  other  mechanisms  (e.g.,  correlated  electrons)  giving  rise  to  the  transport  and 
carrier  concentrations  in  these  compounds  are  presently  ongoing,  and  some  preliminary 
data  was  presented  at  the  recent  MURI  meeting  and  discussed  with  collaborators.  The 
next  steps  include  trying  to  develop  ways  of  controlling  these  properties. 
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2013 


G.  Chen,  “Thermoelectric  Transport  in  Bulk  and  Nanostructured  Materials,”  invited  talk 
at  Gordon  Research  Conference  on  Nanomaterials  for  Applications  in  Energy 
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alloys  and  discovery  of  potassium  as  a  resonant  level”,  International  Conference  on 
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G.J.  Lehr  and  D.T.  Morelli,  “Thermoelectric  Properties  of  Ybi-x(Sc,Er,Lu)xAl3  Solid 
Solutions,”  Poster  presentation,  31st  International  Conference  on  Thermoelectrics, 
Aalborg,  Denmark,  2012. 

G.J.  Lehr  and  D.T.  Morelli,  “A  study  of  the  relationship  between  intermediate  valence 
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Materials  Research  Society,  Boston,  MA,  2012. 

B.  Liao,  K.  Esfarjani  and  G.  Chen,  “First-principles  study  of  thermal  transport  in  FeSb2”, 
ASME  International  Mechanical  Engineering  Congress  and  Exposition,  Houston,  TX, 
Nov.  9-15,2012. 

B.  Liao,  M.  Zebarjadi,  K.  Esfarjani  and  G.  Chen,  “Cloaking  core-shell  nanoparticles  from 
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Thermoelectric  Symposium. 
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G.  Chen,  “Thermal  Transport  and  Properties  in  Nano  structured  Materials,”  Plenary 
lecture  at  18th  Symposium  on  Thermophysical  Properties,  Boulder,  Colorado,  June  24-29, 
2012. 

G.  Chen,  “Nanostructured  Materials  for  Thermoelectric  Energy  Conversion,”  invited  talk, 
ACS  National  Meeting,  Philadelphia,  August  18-23,  2012. 

G.  Chen,  “Nanostructured  Materials  for  Thermoelectric  Energy  Conversion,”  invited  talk, 
Orcas  2012,  International  Conference  on  Energy  Conversion  and  Storage,  Friday  Harbor, 
WA,  September  4-6,  2012. 

Z.F.  Ren,  “Solar  Energy  Conversion  into  Electricity  and  Hot  Water  by  Thermoelectric 
Effect”,  ACP  2012,  Guangzhou,  China,  Nov  7-10,  2012. 
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Z.F.  Ren,  “Aligned  Carbon  Nanotubes  and  Nano  structured  Thermoelectric  Materials”, 
East  Lake  International  Forum  on  Frontiers  of  Science  and  Technology  for  Outstanding 
Overseas  Young  Scholars,  Wuhan,  China,  Oct.  6-8,  2012. 

Z.F.  Ren,  “Nanomaterials,  Physics,  and  Applications  in  Energy  Research”,  Workshop  on 
Energy  and  Advanced  Information  Materials,  Shengzhong  Lake,  Sichuan,  China,  July  22- 
29,  2012. 

Z.F.  Ren,  “Nanostructured  Bulk  High  Temperature  Thermoelectric  Materials”,  MRS 
meeting,  San  Francisco,  CA,  USA,  April  9-13,  2012. 

Z.F.  Ren,  “Electron  and  Phonon  Engineering  in  Nanostructured  Thermoelectric 
Materials”,  The  USC  Smart  Energy  Summit,  Los  Angeles,  CA,  USA,  Jan  27,  2012. 
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H.  Jin,  B.  Wiendlocha,  J.P.  Heremans,  “Indium  is  an  Acceptor  in  Bismuth”,  American 
Physical  Society  March  meeting,  Dallas,  Mar.  2011. 

S.R.  Boona  and  D.T.  Morelli,  “Thermoelectric  properties  of  Cei-xScxPd3,”  Poster 
presentation,  30th  International  Conference  on  Thermoelectrics,  Traverse  City,  MI,  2011 

Y.  Tang,  Poster  presentation,  30th  International  Conference  on  Thermoelectrics,  Traverse 
City,  MI,  2011. 

N.  Heinz,  Poster  presentation,  30th  International  Conference  on  Thermoelectrics,  Traverse 
City,  MI,  2011. 

Z. F.  Ren,  “Nanostructured  Thermoelectric  Materials”,  HI  TEMP  2011  Conference, 
Boston,  Massachusetts,  USA,  Sept  20-22,  2011. 

Z.F.  Ren,  “2nd  Status  meeting  SPP  1386  for  Nanostructured  Thermoelectrics”, 

Wittenburg,  Germany,  July  4-6,  2011. 

Z.F.  Ren,  “Nanostructured  Thermoelectric  Materials  and  Their  Applications  in  Energy”, 
MRS  meeting,  San  Francisco,  USA,  April  25-28,  2011. 
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Z.F.  Ren,  “Nanostructure  Approach  for  High-performance  Thermoelectrics, 
Photovoltaics,  and  Biosensing”,  2011  NES/APS-AAPT  Joint  Spring  Meeting,  UMass 
Lowell,  Massachusetts,  USA,  April  8-9,  2011. 

Z.F.  Ren,  “Nanostructured  Thermoelectric  Materials  and  Their  Applications  in  Energy 
Conversion”,  Symposium  on  Nanophotonics  and  Renewable  Energy,  Beijing,  China,  Jan 
16-18,2011. 

2010 


H.  Jin,  C.M.  Jaworski,  J.P.  Heremans,  “p-type  Bii-xSbx  above  300K”,  American  Physical 
Society  March  meeting,  Portland,  Mar.  2010. 


Colloquium  and  seminar  presentations 

Z.F.  Ren,  Boston  University,  Boston,  Massachusetts,  USA,  “Why  are  nano  materials 
useful  for  enhancing  energy  conversion?”,  November  18,  2011. 

Z.F.  Ren,  University  of  Houston,  Texas,  USA,  “Why  are  nanomaterials  important  to 
energy  physics  and  engineering?”,  November  17,  2011. 

Z.F.  Ren,  University  of  Maryland,  Maryland,  USA,  “Why  are  nano  materials  useful  for 
enhancing  energy  conversion?”,  October  26,  2011. 

Z.F.  Ren,  University  of  North  Carolina  at  Chapel  Hill,  North  Carolina,  USA,  “Why 
nanostructure  is  important  for  future’s  energy  need?”,  September  12,  2011. 

Z.F.  Ren,  University  of  Science  and  Technology  of  China,  Chengdu,  China, 
“Nanostructure  Approach  for  Enhancing  the  Thermoelectric  Performance  of  Bulk 
Materials”,  July  7,  2011. 

Z.F.  Ren,  Minzu  University  of  China,  Beijing,  China,  “Why  nanostructure  is  important 
for  the  future’s  energy  need?”,  June  1,  2011. 


(7)  Plans  for  the  option  years 

So  far,  the  MURI  has  set  the  record  values  of  zT  in  the  cryogenic  temperature  range, 
in  both  n-type  materials  (K-doped  BiSb  alloys)  and  p-type  (CePd2Pt)  materials.  These 
values  will  be  consolidated  (i.e.  assessed  for  durability  and  repeatability),  and  enhanced. 
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Indeed  the  scientific  directions  taken  have  proven  themselves  fruitful,  and  yet  new  ideas 
have  arisen.  In  the  last  two  years,  the  following  research  plans  will  be  enhanced  or 
reduced. 

•  The  work  on  tetradymites  at  Princeton  will  wind  down;  it  will  be  replaced  by 
additional  work  on  the  possible  semiconducting  heavy  Fermion  systems  (FeSb2- 
like)  such  as  RuSb2. 

•  The  work  on  natural  nano  structuring  of  Bii-xSbx  alloys  at  Caltech  will  decrease;  it 
will  be  replaced  by  new  work  developing  a  new  concept  for  cryogenic  cooling 
based  on  the  Thomson  effect. 

•  The  work  on  bulk  Bii-xSbx  alloys  at  OSU  and  MIT  will  be  enhanced,  with  the 
following  focus  points: 

(a)  The  experimental  investigation  of  quantum  structures  (MIT)  has  just  started 
and  will  become  a  focal  point. 

(b)  The  exact  compositions  x  and  K-donor  concentrations  needed  to  develop 
resonances  will  need  to  be  determined. 

(c)  The  effect  of  a  moderate  external  magnetic  field  will  be  investigated. 

(d)  Other  possible  resonant  impurities  will  be  explored. 

(e)  A  new  idea  based  on  the  addition  of  microscopic  (but  not  nanoscale) 
inclusions  of  insulting  materials  will  be  explored,  in  the  hope  of  attempting  to 
prepare  material  with  (1)  Dirac  states  following  the  MIT  calculation  (Dresselhaus- 
group)  (2)  the  idea  of  “cloaked”  core/shell  inclusions  following  the  MIT 
calculations  (Chen-group). 


(8)  Transition  planning 

Given  the  progress  to  date  of  the  values  for  ZT  achieved  by  MURI  team 
members,  contacts  are  being  initiated  with  AFRL  scientists  at  the  Space  Vehicles 
Directorate,  Kirtland  AFB.  Ms.  Gloria  Lehr,  who  is  pursuing  her  Ph.D.  research  at 
Michigan  State  University,  under  the  direction  of  Prof.  Donald  T.  Morelli  under  this 
MURI  program,  and  is  the  author  of  the  work  on  n-type  heavy  Fermion  systems  based  on 
AhYb,  applied  to  the  Air  Force  Research  Laboratory's  Space  Scholar  Program. 

Another  change  being  planned  is  the  move  of  Prof.  Zhifeng  Ren  from  Boston  College  to 
the  University  of  Houston. 


(9)  Issues  of  concern 
None 


(10)  Recommendation  for  funding 

It  is  strongly  recommended  that  this  program  be  fully  funded  for  the  final  2  years. 
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